Breathing Modes and Hidden Symmetry of Trapped Atoms in 2D 
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Atoms confined in a harmonic potential show universal oscillations in 2D. We point out the connec- 
tion of these "breathing" modes to the presence of a hidden symmetry. The underlying symmetry 
SO(2,l), i.e. the two dimensional Lorentz group, allows pulsating solutions to be constructed for 
the interacting quantum system and for the corresponding nonlinear Gross-Pitaevskii equation. We 
point out how this symmetry can be used as a probe for recently proposed experiments of trapped 
atoms in 2D. 
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The problem of Bose-Einstein condensation in an ex- 
ternal potential has received a lot of interest after the 
experimental observation of the condensation in alkali- 
atom vapors JlJ . In real experiments the trapping poten- 
tial is approximately harmonic with frequency u)q. This 
implies specific peculiarities in the behavior of the sys- 
tem. In recent papers a few of these properties have 
been demonstrated. It has been shownjpj that the non- 
linear Gross-Pitaevskii (GP) equation ( |l2| ) for a trapped 
2D system possesses "breathing" oscillatory modes with 
the universal frequency 2loq, describing a pulsation of the 
condensate. The same modes show up in where the 
authors were able to construct explicitly the time evo- 
lution of the GP equation in a time dependent external 
potential in 2D. In || the authors discovered that the en- 
ergy spectrum of a system of trapped particles interacting 
with a 1/r 2 potential is divided in sets of equidistant lev- 
els with the separation 2wo again. To our understanding 
this interesting property has not yet been properly ex- 
plained and the connection between these different sys- 
tems has not been established. 

In this paper we shall show that in these cases the ex- 
istence of these 2u>q oscillations is ensured by a specific 
symmetry property of the system. Proper use of this 
symmetry leads to a transformation which permits a set 
of breathing mode wave functions to be constructed al- 
gebraically not only for the 1/r 2 problem, but also for a 
local interaction in two dimensions. Actually this sym- 
metry is not only a property of the mean field theory 
as found in HQ but of the full quantum theory. We will 
show, that the oscillating solutions of Q are a continuous 
representation of the underlying group 5*0(2, 1). 

To understand the role of symmetry in the problem it 
is useful to consider as an instructive example a system of 
classical particles moving in a harmonic external poten- 
tial Vp t — ^mu^r 2 and interacting with a potential 
V(rj) with the scaling property: V(Xri) = X n V(ri), 

Let us consider the quantity I = J2i r 1 sucn that d t I — 
2 J2i r iPi/ m - Following the usual derivation of the virial 
theorem in classical mechanics we get: 



= 2T - nV - 2V r 
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(1) 



We now see that for a potential with the scaling exponent 
n = — 2 the right hand side of ([!]) takes the form 2E — 
2muj 2 .I 1 where E is the total energy of the system. In 
this case one gets a closed equation for /: 



dfl = -Au%I + AE/r 



(2) 



with the obvious solution: I = Acos(2uiot + 7) + 
E/(i7iu}q). Thus the existence of the "2u;o" modes is con- 
nected with the n = — 2 scaling of the interaction poten- 
tial. In a 3D system the only potential possessing this 
property is the 1/r 2 interaction used in S. But in the 
quantum 2D case the Fermi "pseudo-potential" \g5 2 {r) 
as used in the GP equation gives the same scaling. 

The equations above can be rewritten by introduc- 
ing (the notations will become obvious later) L + = 
J2iPi r i/ 2 + i( E - mw$I) /(2w ). We simply get d t L+ = 
i2ujQL + . The phase of L + varies linearly with time: 



$(pi,n) = — — ImlogL + 

ZUJQ 
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■ arctan 



(E - mLU 2 I)/(2io ) 



2^o V, /',''- 2 

$(Pi,ri)-$(p?,r?)=t-to, 

where P°,t° are the coordinates at t = to. $ can now 
be used to determine the "abbreviated action" S(E,ri) 
which is a function of the energy and the coordinates at 
the end of a path. It is determined by the Hamilton- 
Jacobi equation H(d ri S,ri) = E with t = &eS = 
3>(d ri S,r). In hyper-spherical coordinates in the space 
of all Ti with r = vJ, $ is only a function of r and d r S. 
Therefore we have S(E, ri) — S(E, r) + So (a*) where the 
Q.i denote all other coordinates, clearly showing that the 
coordinate r ox I totally separates. 

One of the most powerful methods in physics is the use 
of symmetries and groups. One way is to use the invari- 
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ance of the Hamilton or the action under certain trans- 
formations, another is recognizing that the Hamiltonian 
is a part of some larger algebra. The most famous text- 
book example is the algebraic solution of the harmonic 
oscillator using the spectrum generating Hciscnberg al- 
gebra [H, a^] = ±ujQa^ (we put % = 1 throughout the 
paper). 

We will now discuss such a spectrum generating sym- 
metry for the (now quantum-mechanical) problem of in- 
teracting particles in a harmonic trap. First we will con- 
sider the effect of a scaling-transformations for the Hamil- 
tonian H = ~ 5 kA l + Ei<j v ( r ^ 
external potential: 



r,) without an 



r -► Ar,*(r) -► X d/2 ^(\r) 



(3) 
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Ho is scale invariant if V(Xr) = V(r)/X 2 . This is the 
case for an interaction of the form V(r) — g/r 2 in any 
dimension but as mentioned above also for 



V(r 



\ 9 5 2 {r-r') 



(4) 



in two dimensions. The first case is known as the 
Calogero-Sutherland model § and widely investigated. 
All our results apply for both cases as we will only use 
symmetry properties connected to scale transformations, 
but we concentrate on the problem of bosons with a local 
interaction and only shortly comment on the relation to 
known results for the Calogero-Sutherland model. Actu- 
ally the local interaction is the most important case, at 
least for neutral atoms at low energies, where the range 
of interaction is small compared to all other scales. The 
statistics of the particles does not play a role. 

It is important to note that a (5-function interaction is 
not well defined in two dimensions due to logarithmic ul- 
traviolet divergences which are cut off by the finite range 
ao of the interaction. This length oo obviously breaks 
scale invariance of Ho and will therefore modify our re- 
sults. Nevertheless this effect will be small as long as 
ao is smaller than any other scale in the system. This 
is most clearly seen in the classical wave-limit (i.e. the 
GP-equation) which we will discuss later, where such a 
problem is absent. 

Adding an external potential H = H a + H pot , H pot = 
J2i \ mu >o r 'i obviously breaks scale invariance, as i? po t — > 
X 2 H pot under a scale transformation. However, due to a 
special property of the harmonic oscillator, still a power- 
ful spectrum generating symmetry exists. 

The important step is to recognize that the commuta- 
tor of the harmonic potential with the Hamiltonian or the 
time derivative of J>2% r 1 lSi proportional to the generator 
of scale transformations. 

lH po uH} = lY / 1 r ^ 2 r 2 ,Y,-^] 



= ^2^l(V l r l + r l V i ) =iulQ (5) 

i 

Q = Y1 ^(Pin+riPi). 

i 

Q is the generator of scale transformations as it describes 
the translation of the coordinates by an amount pro- 
portional to 7"j. 

We can collect our results (0,0) in the following algebra: 



[Q,H ] = 2iH , [Q,H pot ] = -2iH pot , [H pot ,H] = iu 2 Q 



or using 

Li 

we get the algebra: 
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— { H + H pot ) = —H 



(6) 



[Li,L 2 ] = -iL 3 , [L 2 , L 3 ] = iLi, [L 3 , L x ] = iL 2 . (7) 

This is the well known algebra of SU(l, 1) or SO(2, 1), 
the two dimensional Lorentz group. L\ and L 2 are the 
generators of the two "boosts" and L3 = ^j-H, i.e. the 
generator of time-translations, is the analog of the gen- 
erator of the rotation. With i ± = -^(Li ± iL 2 ) this 
reads 

[H,L ± }=±2u L ± ,[L+,L-} = -^-H (8) 

Note the minus-sign in the last equation, indicating that 
the group is the Lorentz group SO(2, 1) (or SU(1, 1)) and 
not 50(3). 

One important consequence of this spectrum gener- 
ating symmetry is that the Hilbert-space will separate 
into irreducible representations of the group. If the en- 
ergy is bounded from below these are discrete infinite- 
dimensional representations with no upper bound. Start- 
ing from the lowest eigenstate in one of the representa- 
tions with energy E , H\^ ) = ^bl^o) one can construct 
higher states with energies E + n2u} , n — 1,2,... by 
applying L+ (use HL+\^ ) = (L+ H + 2lu L+)\^ ) = 
(Eo+2lo )L + \^ )). |*o) is annihilated by L~ . Obviously 
an infinite number of excitations with energies n2wo ex- 
ists which we will identify with the breathing modes of 
the system. 

Also the time dependence of all the operators, which 
are part of the algebra can be given explicitly: 



(9) 



Defining the operator for the mean square displacement 
of the particles I(t) — J2 r i (f) one finds 
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(10) 
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This equation or the corresponding differential equation 
dfl = — (2uj ) 2 (I — H/(muj 2 )), coinciding with the classi- 
cal Eq. (||) , clearly show that the variable I is separated 
from all other variables, in analogy to the center of mass 
motion. This separation is known in the case of the 1/r 2 
interaction As the frequency lu of the external po- 
tential obviously couples only to this coordinate also the 
full nonlinear response to a change of the potential can 
be calculated ||,f§]. 

For the expectation value we directly get 



I(t) = (I(t)) =I + Acos(2u t + 7) 



(11) 



with I Q = (H)/(mu$) and Ae~^ = % /2/(mcj )(i + )t=o- 
The same solution was found recently for the GP equa- 
tion in D = 2 HQ], but here we show that it also holds 
for the full quantum system. It is valid not only for the 
expectation value but directly for the operators (|l^) and 
is due to a simple underlying symmetry. We will discuss 
the change of this equation under a group transformation 
later. 

The Casimir operator of this group C = H 2 — 
(2u>o) 2 (L 2 + L 2 ) commutes not only with the Hamilto- 
nian, giving a new conserved quantity, but also with the 
generators of the group and is the analog of the invariant 
line element in special relativity. 

The classical wave limit of interacting bosons, i.e. the 
non-linear GP equation ||, has recently attracted con- 
siderable attention as it accurately describes the Bose- 
condensate of trapped atoms and allows for reliable an- 
alytic and especially numerical calculations Q. It is 
also an interesting problem on its own for mathemati- 
cal physics. We will show that the previously discovered 
analytic solutions Q actually form a continuous repre- 
sentation of the group SO(2, 1) for the GP-equation. 

In Q the following transformation of the solution of 
the GP equation (m = h = 1) 



= -iA* + ^ 2 r 2 *+.g|*| 2 *-/i*. 



was considered (in our notation): 



V 2 (r,t) = exp[i(b + cr 2 )} — *i(m,t), 



u = 



-,r(t) = 



df 



(12) 

(13) 
(14) 



If the functions a(t), b(t) and c(t) fulfill (at — dta etc.): 



a, 



<iac, bt = I 1 

a 



a u a - — 



(2Q ) 2 



(2w ) 



= 0, 



(15) 



and if *f?i is a solution of ( |l2] ) then \&2 is also a solution 
of the GP-equation with a possibly time dependent fre- 



quency Qo(t). Differentiating ( |l5| ) gives the linear equa- 
tion aut + (2tio) 2 a>t = — 4(2>o^ot a which once again demon- 
strates the universal nature of these modes; the initial 
values must fulfill (|l~5|). 

We will now consider the case ujq — ^0 — const where 
the differential equations ( |TJ ) can can be solved directly: 



a = sinh r\ cos(2o;ot + 7) + cosh r\. 
t= — (arctan (e~'' tan (uiot)) + irn) 
1 



(16) 



/ sinh 77 + cosh n cos 2u> t \ 

. arccos — sign(sm 2ujot) 

2uj ^ V cosn V + smn V cos 2^0* / 



27rn 



(17) 



and c(ji,t),b(ji,t) accordingly. The integer n has to be 
chosen to get a continuous solution, n(t) = [u t/n+ 1/2], 
where [x] denotes the greatest integer less than x. A par- 
ticularly important example is the case when the initial 
solution is a static one, e.g. the groundstate, = 1 i'a{ r )- 
The transformation builds then from such a static solu- 
tion an oscillating breathing solution. In this case the 
parameter 77 of the transformation defines the relation of 
the energy of the new solution E to the static one Eq 
according to E = cosh(rf)Ea- In the following we set the 
phase 7 = in (|l^) as a finite 7 can be achieved by a 
simple translation in time. 

Together with the time translations the discussed 
transformation form a group, indeed a continuous rep- 
resentation of 5*0(2, 1). A general group transformation 
J7t2,n,ti can be described by a initial translation in time 
by fx, the above described scaling transformation param- 
eterized by 7/ and a second final translation backwards in 
time by ti. This is the analog to the description of rota- 
tions by Euler angles. With r' = r(r\, t — t\) and a', b' , d 
accordingly we have 

EWi* = e« b ' +c ' r2 ^y(-^=r,T> + t 2 ). (18) 



It is now simple to work out the multiplication rules 
of this group - they are the same as for two dimensional 
Lorentz transformations. For example if one performs 
two successive "Lorentz boosts" in the same direction 
one has to add the rapidities: 



Ui 



0,772,0 ^0,771,0 



= u. 



0,7ji+r7 2 ,0- 



(19) 



We can see this by calculating the corresponding function 
a(t) for two transformations using (|l7|): 

a — (cosh 772 + sinh 772 cos 2ujot) 

x (cosh?7i + sinh 771 cos(2a>oT(>72) t))) 
= cosh(77x + 772) + sinh(7?i + 772) cos 2tu t. (20) 

To make the connection to the previously constructed 
algebra, one has to calculate the generators of the group, 
i.e. the infinitesimal transformations 5rj — > 0. 
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1 + St] cos 2cjgi, T 



t — Sn sin 2u>ot- 

2iv Q 



Under such a transformation the wave function changes 
by 5V = C/ ,5,,,o* - 



<5* 



1, 

- A {rp 



■ pr) cos 2wo^ 



1 • „ ( 9 
- — sin2w t i — 

2w V ^ 



2-w, 



2„2 



* (21) 



with p 



= — iV. This has to be compared with ii = 
pr) and L 2 = 2^(# - 2i/ pot ). Identifying 
with idt and noting that in the Heisenberg picture 
L 2 (t) = V2ImL+(t) = L 2 cos 2tJoi — Li sin 2wo^ we can 
identify £/o,?7,o with the transformation generated by 
Accordingly L\ generates the transformation f/o,?7,7r/(2wo) 
and iJ obviously J7t,o,o- Actually the solutions generated 
by Li and L 2 describe the breathing of the system, i.e. 
a pulsating motion. 

To finish this section we directly evaluate the effect of 
a group transformation to the mean square displacement 



of the bosons I(t) = J r 2 \^\d 2 r. From (jl7|) we know, 
that for any given solution \l/o, lit) is of the form I{t) = 
C(coshr/o + sinh?7o cos2woi) as long as the system does 
not collapse, i.e. as long as I(t) > 0. (Such a collapse is 
of course possible only for an attractive interaction, i. e. 
for g < 0.) With = C/ ,,,,o* and using @) we get: 



I'(t)= J r 2 \*'(r,t)\d 2 r = a(t)I(T(t)) 



= C(cosh(?7o + if) + sinh(?7o + 77) cos 2w i). (22) 

In particular, any solution without a breathing motion 
(r/o = 0) will be transformed to a breathing one and on 
the other side any breathing motion can be transformed 
to by choosing 77 = —770. 

As mentioned below Eq. ( pi] ) all these properties are 
a direct consequence of the underlying group structure 
and aj 2 Ccosh(?7o + i]) can be identified with the en- 
ergy of the system, while the amplitude of the oscilla- 
tions C sinh(?7o + 77) is 2/w -\/ (£i) 2 + (L 2 ). The quantity 
C = (H) 2 - (2uj q ) 2 ((L 1 ) 2 + (L 2 ) 2 ) is conserved under a 
group-transformation (cosh 2 — sinh 2 = 1) and plays the 
role of the Casimir invariant. Generally, if C is posi- 
tive, one can always find a group transformation, so that 
(L[) 2 = (L' 2 ) 2 = and therefore I(t)' = = 
const. For C < this is not possible and I(t) describes 
the collapse of the system as has been discussed in M . 

To conclude, we have shown that the existence of 
breathing oscillations of 2D atoms in a harmonic trap 
is ensured by a hidden symmetry of the system. If 
the atoms interact to a good approximation by a Fermi 
pseudo-potential, i.e. by a local interaction, we expect 
well defined modes with a frequency of exactly 2tu . We 
identified the corresponding solutions of Q with contin- 
uous representation of the underlying "Lorentz" group 
SO(2,l). 



The experimental conditions for the realization of a 
2D system in a magnetic trap look quite promising in 
the case when most of the atoms are in the condensate. 
The system can be considered as 2D if lo z ojq is ful- 
filled. But our calculations are also valid for a system 
without condensation, e.g. for fermions. In this case the 
condition for two-dimensionality also demands that the 
effective temperature is less than the level separation in 
z-direction: T huj z . 

On contrary one can also use a in z-direction prolonged 
trap H and excite oscillations in the condensate with no 
z-dependence. 

Important new possibilities are opened by a recent pro- 
posal ||] to confine atoms in a 2D optical dipole trap. In 
such an experiment the appearance of sharp 2wo frequen- 
cies would be the first demonstration of the two dimen- 
sional nature of the system. Experimentally it is also easy 
to excite the 2ujq- modes by a change of the external po- 
tential as described in Q| . Precision measurements of this 
response permit to check the assumptions of the model. 
Especially the validity of a local interaction, which was 
generally accepted so far, could be investigated. It will 
be interesting to study deviations from such a description 
both experimentally and theoretically. 
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